ON MATRIX SCHRODINGER UNITARY GROUPS IN 
PARTICULAR REPRESENTATIONS OF FINITE DIMENSIONAL 
QUANTUM DYNAMICAL SYSTEMS 

^ ! 
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' Abstract. In this paper we study some particular types of matrix Schrodinger 

unitary groups of the form exp(— itH) where H S Mjv(C) is the Hamiltonian of 
a given quantum dynamical system modeled in the finite dimensional Hilbert 
^\ ■ space "H. Once we have defined a particular matrix Schrodinger unitary group 

we perform some estimates for its approximation and its corresponding im- 
plementation in the numerical solution of the finite dimensional Schrodinger 
' evolution equation to that it is related. 

■ 1. Introduction 

In this work we will focus our attention in the study some properties of Ma- 
trix Schrodinger semigroups that will be described in general by the set {Si :— 
t:^ ' exp(— iffl) : t £ R}, where H e ^{Ti.) is the hamiltonian of a prescribed quantum 

^ . dynamical system and T-L := 'H{G) is a finite dimensional Hilbert space related in 

some suitable sense to a weighted graph G := {Vg^Eq) that will be called repre- 
sentative graph, the elements of the semigroup clearly satisfy the conditions: (i) 
; Sq = 1«, (ii) St • S,(-) = St+,(-) and (iii) lim;,^o+ ^h<i> = 0,V0 G I? C be- 

• . sides the condtion (iv) ||St'!/'o|| = IIV'oH, will be also satisfied when H is self adjoint. 

' The Hamiltonian H e J/f^H) presented above is related to a prescribed quantum 

dynamical system trough the evolution equation given by: 



X 



Em)) = '^\m) (1-1) 

with ip{0) — ipQ E H and where E — > i/HDt, here H e ^{H) will in general have 
the form H = Hq + Hi with Mq € Mn{C) self-adjoint, and with Hi e Mw(C) 
^ . diagonal, for simplicity, in this work we will consider our scale such that h = 1. 

In the following sections we will implement some operator theory techniques in 
the theoretical analysis of the approximation schemes for the matrix Schrodinger 
semigroups and in the end some numerical implementations will be presented. 

2. Particular Representations of Finite Dimensional Quantum 

Dynamical Systems 

In this section we will present the main ideas behind particular representation 
techniques of finite dymensional quantum dynamical systems. 
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2.1. Generalized Matrix Form of Kets and Bras. In this work wc will consider 
that the space of states H oi a. finite dimensional quantum dynamical system, is a 
finite dimensional Hilber space, whose inner product is induced by an inner product 
matrix Mh S Mjv(C), wich is a symmetric positive definite matrix that satisfies 
the relation 

{u\v) := ti*Muv (2.1) 
since Mh is symmetric positive definite, we can obtain a factorization of the form 
Mh := W^Wh, where Wh e Mjv(C) represents the formal square root of Mh- 
From the above relations we can obtain matrix representations for kets, bras and 
induced norm by (•!•) operations according to the rules: 

\u) i — > u (2.2) 

{u\ ^ u*Mh (2.3) 

^ (u|«)^/^ (2.4) 

from the definition of the norm operation, and if we denote by (•, •)2 the usual inner 
product in given by {x,y)2 ■= y*x it can be seen that {x\y) = {y,x)2 and also 



that 



2 

||u|| := (WhW, Wmu)1^^ = ||Wh«||2 , (2.5) 
in a similar way the induced matrix norm in Mjv(C) by ||-|| can be expressed in the 
form 

||A|| := sup ||Aw|| = IIWaAWii^ll^; (2.6) 

||u|| = l 

from the relation between (•!•) and (•,-)2) it can be seen that for a given matrix 
A e Mjv(C), one can compute its adjoint At with respect to (-I-) using the following 
expression 

A'f := M-^A*Mm. (2.7) 

2.2. Particular Ladder Operators. For a given finite dimensional quantum dy- 
namical system with space of states H := ?^(G), one can find or obtain an or- 
thonormal system X := {\k) , 1 < fc < N} C H with respect to (-I-) that will be 
called particular analysis basis, and that is related to a given representative graph 
G := {Vg,Eg) with Vq := {I,-" i-^} and Eq '■= {{{ij) according to the 

rule 

k^\k). (2.8) 

Remark 2.1. It is important to note that an important idea behind the particular 
analysis basis, is to find a basis that provides some advantage for the analysis of 
the quantum dynamical system in wich space of states we define it. 

Once we get an orthonormal system X C "H it is possible to obtain a matrix 
N e Mjv(C) that satisfies N\n) = n\n) and will be defined by 

N 

N:=J2k\k){k\ (2.9) 
fe=i 

in a similar way one can define two particular ladder operators at, a € Mjv(C) 
defined implicitly by 



at|n) := ^/[n + l\N]\[n + 1\N]) (2.10) 
a\n) := ^/[^] \[n - 1\N]) (2.11) 
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where [p\q\ := 1 + {{p — 1) mod q). From the imphcit definition of the particular 
Ladder operators we can obtain the following exphcit definitions 



a 



N-l 

t ■= |1)(7V|+ ^ VfcTT|fc+l)(A:| (2.12) 

k=l 
N 

a := |7V) (1| +^ Vfc|fc- 1) (/fc| . (2.13) 

k=2 

It can be seen that N = a^a. Since we will have that X C H is an orthonormal 
system, also that hJ = Hq, and \m) = (TO!)-i/2(a1')™ |1), and if we take Ei := 
|1) (1|, then we can express Hp £ Mn{C) in the form: 

Ho := ^(fc|Ho|/)|fc)(/| (2.14) 

(ki) 

= ^ (M^(«t)'ciE^„fc + ^ »2|) ((„t)fe]E^„^ + (at)'Eia^-l2.16) 

(fc) {k>i) 

that will be called particular representation of Hq G Mjv(C) with respect to X C 



3. Time Evolution and Matrix Schrodinger Unitary Groups 

For a finite dimensional quantum dynamical system with state of spaces 'H(G), 
whose time evolution is modeled by the equation 

iW{t))=mm)) .o1^ 

\m) = iV'o) ^ ' 

one can obtain an expression for its wave function \ip{t)) using the matrix valued 
function C/f : R — >'H:t\ — > e"*", in the form 

|V(i)> e-^« IV-o) (3.2) 

in many applications H G A/jv(C) is a structured matrix obtained using several 
matrix operations between matrices of lower order. In the following subsection 
we will consider particularly important cases of interacting and non-interacting 
quantum systems. 



3.1. Matrix Hamiltonians. For a sequence of finite dimesional Hilbert spa-ces 
of the form {'Ha]a=i: if '^^'^ obtain a particular analysis bases sequence {^a}, 
such that the non-interacting Hamiltonian of each finite quantum dynamical system 
whose space of states is defined by Ha can be particularly represented by Hq, and 
defined according to (|2.14l) . obtaining a sequence {Hq} that can be used to compute 
the non- interacting part an interacting Hamiltonian H g Mjv(C) using the following 
expression 

Ho :=0Hq = ^{H„}0^" (3.3) 

(a) (a) 
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where {cq := {Sr,a)r}a C (Z(j")^^ is the canonical basis for the space of muhi- 
indexes of length M £ Zq, and where in general for any multiindex r e (Z^j")^^, 
{H^}®'' is defined by 

{Ba}®'':=(g)B;= (3.4) 

(a) 

here Iq denotes the identity matrix in "Hq. When we want to compute a solution 
to (|3.ip and if H := Ho + Hi, where Hp is an non-interacting hamiltonian and Hi 
is diagonal, it is useful to compute an integrating factor that is a solution to the 
matrix differenial equation 

Uo = 1 ^^-^^ 

A solution to this equation will have the form Ut ■= e"*™", and clearly the set 
{Ut,t € R} will be an unitary group of operators, now this integrating factor can 
be used to solve (|3.ip using the formula 

m)) = Ut \^o) + f dsUt-s^m)) (3.6) 
Jo 

whose solvability have been discussed in 9i, in particular the case where both the 
interacting and non-interacting parts of the hamiltonian H can be expressed in a 
non-interacting Hamiltonians form, i.e., H = Mq, it is important because we can 
express the solution to p.ip in the form 

\m) Ut iV^o) = e-™" iVo) (3.7) 

where 

\m) <S)M)) (3-8) 

(«> 

IV'o) := (g)|V^o,a). (3.9) 

(a) 

It is also important to note that for some complex systems modeled by non- 
interacting type Hamiltonians, even if their wave funtion can be expressed like 
([XT]) , the computation of e"*"» that can be expressed in the form 

e-™° := (g) e-™° (3.10) 

can become a hard computational problem, this is the reason to implement some 
numerical techniques that will be presented in the next section. 

3.2. Approximation of Matrix Schrodinger Unitary Groups. When we want 
to compute an approximation of a particular matrix Schrodinger unitary group 
{Ut '■— e^'*'^,t e R}, we can start approximating Ut in [0,t] C R, for a given 
1 > T := ft,t/||Hj| e R, with H 3 ht < I, this restriction for t will help to ensure 
that the sum 

k=0 
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remains bounded with respect to |j-||. Clearly this sum represents the first m terms 
of the Taylor polynomial of Ut around t — 0, now if we take the Pade representation 
of the approximant U we obtain the following expression: 



with 



{p + q- 3)¥- 



Npqi-ITM) 

Dpgi-trtt) 



E 

i=o 

E 

3=0 



{p + qy.jl{p - j)l 

{p + qWiq- j)\ 



(3.12) 

(3.13) 
(3.14) 



It can be seen that taking S :— Npp{—iTM)^ we will have that Dpp(— irH) = S^', 
and if we take §+ := (S^^)^^, then we can express p.l2p in the form 



U = §+§ 



(3.15) 



From the relation of p.l2p with the Taylor expansion of e and the Picard's 
resctriction for r £ R+ in (|3.1ip . we can obtain the following estimate 



Lemma 3.1. ||e 

Proof. 



(2p+l)! C,p,2p+1 



E (;^-Cp,fc)(-*rH)' 



< 



< 



< 



fc=2p+l 

oo 

E 

fc=2p+l 
1 



fc! "'^'^ 



(2p+l)! 
1 

(2p+l)! " 



r2P+l llml|2P+l 



(3.16) 

(3.17) 

(3.18) 
(3.19) 

□ 



Since H will be considered in general self adjoint, i.e., ~ H, we will have that 
HI is normal, hence can be factored in the form H = VBV*, with D := diag{c?i}, 
and taking A := i?pp(— zrB) we obtain 



U = §+§ 

= VA+V*VAV* 
= VA+AV* 



(3.20) 
(3.21) 
(3.22) 



wich implies the following result. 



Lemma 3.2. U*U = 1. 
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Proof. 



U*U = VA*A"iV*VA+AV* 
= VA*A-iA+AV* 
= VA*A+A-iAV* 



From the above relations we can see that the operator 

U := Wg^UWH 

satisfies the following relations 

= ||WHWgiUWH0|L = ||UWh0||2 = ||Wh0||2 = 



(3.23) 
(3.24) 
(3.25) 
(3.26) 
(3.27) 
□ 

(3.28) 
(3.29) 



wich implies that ||1LJ|| = 1, the adjoint of U can be obtained in the following way: 

= Mh^1LJ*Mh (3.30) 
= M-1(W-1UWh)*Mh (3.31) 
;iW^U*W+M„ (3.32) 

(3.33) 
(3.34) 

the matrix Schrodinger unitary group relative to H, will have the form {Ufc := 
{]''{■), k e Z+}. It can be seen that 

U^U = 1 = UU^ (3.35) 

and this implies that the discrete time reversal Schrodinger unitary group will be 
given by {tS-k ■= k G Z+} and will be coherent with the local time 

reversibility of Schrodinger unitary groups. From lemma L 13. H and taking the time 
interval [0,r7ir] C R we can obtain the following. 



Lemma 3.3. 



-zmrM 



- (2p+l)!"'r 



Proof. If we denote by Cp^k the multinomial Fade coefficients, then we will get 



-imrM 



U" 



(— ztotI 



< 



< 



< 



fc! 

,771*-' 



fc=0 



E 

oo 



fc=2p+l 

(2p+l)! - ■ 



(3.36) 
(3.37) 

(3.38) 

(3.39) 

□ 



MATRIX SCHRODINGER UNITARY GROUPS 



7 



Now if we denote by {Uq} the sequence of matrix approximants corresponding 
to the sequence of matrix Schrodinger basics {e"*"^^"}, and if we take {ha} to 
be the sequence of basic time step sizes considered for the approximation of each 
matrix Schrodinger basics, and if we take ^-'^"^'^^o defined according to (I3.10p and 
the non-interacting approximant U € Mn{C) defined by 



U:=(g)U„ 

we can obtain the following estimate. 

Theorem 3.1. There exists h £ (0, 1) C such that 
l)m^P+^/{2p+iy.h^P+\ 



(3.40) 



g— zmrMo _ "0m 



<(2 



M 



Proof. If wc define h :— sup„ ha and take A^.q := (e^""'^ 
operation 

M 

{A,}0{BJ (g)(A„(ej, • rp) + Mp{ep ■ Sp)) 
p=i 

where r,s £ (Zj)*^ are multi-indexes that that satisfy = | 

M,r-s = and |r| > 0, to express e"™^""" - U"* in the form 



(a) (a) 

^ {A„,j0{ur} 

{\r\ + \s\=M) 



U™) we can use the 



= 1, \s\ 



(3.41) 
(3.42) 



Therefore, 



^ {A„,.}0{ur} 

{\r\ + \s\=M) 
(|r|<M) (r) 



< 



< 



V T\— h^P+^ 

(|r|<Af) (r) ^ ^ ' 

(2^^-l)m2p+i^,p^, 



(2p+l)! 



wich provides the desired result. 



(3.43) 
(3.44) 

(3.45) 

(3.46) 

□ 



4. Numerical Examples 

In this section we present some basic examples to ilustrate the implementation 
of the ideas presented here to compute the evolution of states in a particular finite 
dimensional quantum dynamical system. 

Example 4.1. For a double-slit experiment under simple-absorption photonic con- 
ditions, modeled by a quantum dynamical system with a space of states H ^ , 
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and described by the representative graph G := [VcEg) with vertex and edge sets 
given by 

Vg := {1,2,3,4,5} 

Eg := {((1; 2) , 1), ((1; 3) , 1), ((2; 4) , 1), ((2; 5) , 1), ((3; 4) , 1), ((3; 5) , 1)} 
that is isomorphic to If we take X := {ei~} C "H with Ck ■= {5k,j)j> then we will 



Figure 1. Representative Graph G := {Vg,Eg) of a Double-Slit Experiment. 




have that 
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(4.1) 



(4.2) 



(4.3) 



clearly N = a^a, now, if the Hamiltonian H e M5(C) of the system is defined by 



:A(G) := 



/ 1 1 \ 
10 11 
10 11 
110 
\ 1 1 / 



(4.4) 
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then we will have that e 
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(4.5) 



(4.6) 



also we will have that the Third order Fade approximant that coincides luith the 
Crank-Nicholson scheme will have the form 

/ T^+2 2ri 2 r i 2 2 \ 



u 



V 
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(4.7) 



3t2+2 

by L \3.3\ we will have that 
N G M5(C) con be used to compute the expected state of the system using the 
expression [(N)^,] , with [q] := {p Q Z : \q — p\ ~ mm{\q — r\,r € Vq}} and where 

(N)^ := (V^olUl-NUfe \iPo) = (VolUla^alJfe IV'o) = ||alJfeVo|, 
where \^q) is the initial state of the system. 



(4.8) 



Example 4.2. In this example we will consider a quantum system with three par- 
ticles that evolve in a Fock T-L®'^ space based on the space of states of the above 
example, i.e., T-L®^ := 7i(E)H<E)'H, also we will consider that X®^ -.= {\ijk) ,i,j, k e 
{1, • • • , 5}} and that the hamiltonian of the system will be an interacting hamilton- 
ian, in this particular case will be given by H := H0IHI0H + WqI®"^, where G C 
is an absorption constant related to the media where the system evolves, the group 
basics e"*"^^ and its third order Fade approximant will be given by 

,,3 , _i^ra _ -■.'rH ^ „-iTU\ 



U := {l-iTUjl)-^{l+iTu:l)\]®\]®\] 



(4.9) 
(4.10) 

respectively, this example is not very complex yet, but even this simple example 
provides an idea of how useful can be T \S.1\ to obtain an estimate of the form 



-zmrM 



U" 



< 



-h'^ . In this case we can also find an expression for the 



expected state of the system that will have the form 

LL(N),11 :^ L(Ni>felL(N2)felL(N3).l 

= Lii«wo,iinLii«uv'o,2iiniiauv^o,3in 

where jf/'o) = IV'o.i) ® IV'o, 3) ® IV'o, 3) is the initial state of the system. 
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